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Abstract 

We demonstrate that the Fokker-Planck equation can be generahzed 
into a 'Fractional Fokker-Planck' equation, i.e. an equation which includes 
fractional space differentiations, in order to encompass the wide class of 
anomalous diffusions due to a Levy stable stochastic forcing. A precise 
determination of this equation is obtained by substituting a Levy stable 
source to the classical gaussian one in the Langevin equation. This yields 
not only the anomalous diffusion coefficient, but a non trivial fractional 
operator which corresponds to the possible asymmetry of the Levy sta- 
ble source. Both of them cannot be obtained by scaling arguments. The 
(mono-) scaling behaviors of the Fractional Fokker-Planck equation and 
of its solutions are analysed and a generalization of the Einstein relation 
for the anomalous diffusion coefficient is obtained. 

This generalization yields a straightforward physical interpretation of the 
parameters of Levy stable distributions. Furthermore, with the help of 
important examples, we show the applicability of the Fractional Fokker- 
Planck equation in physics. 
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1 Introduction 



The Fokker-Planck equation is one of the classical, widely used equations of 
statistical physics. It describes a broad spectrum of problems related to the 
evolution of various dynamic systems under the influence of stochastic forces 
and has numerous applications, see, e.g. Usually, the Fokker-Planck equa- 
tion can be derived following the Langevin approach, that is, starting from 
the stochastic " equation of motion" for the dynamic variable whose probability 
distribution we are interested. In this approach, the basic assumptions on the 
"random force/source" in this equation of motion are usually that they have: (i) 
Gaussian statistics, and (ii) delta-correlated correlation. These two assumptions 
yield the (classical) Fokker-Planck equation. 

These assumptions are physically motivated by the fact that the random 
source is a sum of a large number of independent identical random "pulses". If 
these quantities possess a finite variance, then, according to the Central Limit 
Theorem, the distribution of their sum tends to the normal law when the number 
of pulses go to infinity. 

However, the Central Limit Theorem can be generalized for independent 
identically distributed (i.i.d.) random variables having non finite variance. In- 
deed, Levy and Khintchinc ||, ||, ^, |j discovered a broader class of stable dis- 
tributions. They correspond to the limit of normalized sums of i.i.d. stochastic 
variables. Each stable law has a characteristic index a (0 < a < 2), often 
called the Levy stability index or the Levy index, which is the critical order for 
the convergence of statistical moments. Indeed, a statistical moment of a given 
stable law is finite only if its order /i is strictly smaller than its Levy index a (i.e. 
^ < a). Every moment of order higher order (including: fi = a) are infinite or, 
as often said, divergent. The only exception is the normal distribution which 
corresponds to the particular stable law which has its Levy's index a — 2 and 
the exceptional property that all its moments are finite. 

The classical and rather academic example of the application of the Levy 
stable laws is the Holtzmark distribution |^ which is the distribution function 
of the gravitational force created at a randomly chosen point by a given system 
of stars. It is assumed that the system of stars is a (statistically) homogeneous 
set of physical points which mutually interact according to the gravitation law. 
Using dimensional consideration it can be shown Q with a rather straightfor- 
ward scaling argument that this distribution corresponds to a stable law with 
index a = 3/2. We will see (Sect.^ that the following developments allow to 
generalize broadly this result. 

Some other examples of application in physics of stable distributions can be 
found in review papers |^ and in references therein. 

However, let us recall that presumably the most well known application of 
Levy stable laws in physics corresponds to the anomalous diffusion associated 
with a Levy motion, also often called a 'Levy fiight' Indeed, one expects 
that with a Levy stable forcing the cloud of particles will spread much faster 
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(for large times: t » 1) than for a brownian motion. More precisely, we will 
confirm that the radius r{t) of the cloud, at time t, has the following scaling 
law: 

r{t) - t^/" (1) 

the lower bound being reached for the normal diffusion {a — 2), Figs.l and 2 
display illustrations for comparison. 

This scaling relation (Eqj^) is obviously incompatible with the Fokker-Planck 
equation, unless one substitutes a formal (a/2) — th power of the Laplacian to 
the classical Laplacian, therefore considers a Fractional Fokker-Planck equation, 
as suggested by Several authors, following different approaches, considered 
generalization of the Fokker-Planck equation in order to encompass the Levy 
anomalous diffusion. On the one hand, particular cases of the Fractional Fokker- 
Planck equation were obtained ||l2, 13, However, this is not the only way to 



generalize the Fokker-Planck equation in order to respect the anomalous scaling 
relation of (Eq|l[ Indeed, based on Tsallisjl^'s generalization of statistical 
mechanics, a nonlinear Fokker-Planck equation JTs] , [l6| has been introduced 
and it was demonstrated |l^ that its solution respects also Eq.|l|. These 
approaches will be discussed and compared to ours in Sect. |^. 

However, the interest in Levy laws is not limited to anomalous diffusion. For 
instance, the rather large subclass of extremal 1// or "pink" Levy noises have 
been attracting much attention in the framework of multifractal fields. Indeed, 
it was shown that they correspond ||l9| to the attractive generators of " universal 
multifractals" , which are the limit processes, under rather general conditions, 
of nonlinearly interacting i.i.d. multifractal processes. It is worthwhile to note 
that different techniques have been developed to simulate |2^, |2^ or analyse 
[ p2| psi multifractal fields within the framework of Universal multifractals. Let 
us emphasize that the rather straightforward " Double Trace Moment" technique 
yields rather directly an estimate of the Levy index a of the generator. 

Furthermore, in most recent developments of multifractal studies, in par- 
ticular those related to predictability of multifractal processes , one needs 
having a kinetic equation for the generator, because the orientation of time axis 
becomes essential, contrary to earlier simulations, where the generator was ob- 
tained by isotropic (in time as well in space) fractional integration over a white 
noise. 

As there is the need for a kinetic equation in the context of other exam- 
ples/applications of Levy laws in physics, this paper is devoted to establishing 
the corresponding "Fractional Fokker-Planck" equation. 

In this paper we consider the time evolution of a stochastic variable forced 
by a random source having a stable distribution, i.e. a generalized Langevin 



equation (Sect. 2.1). We derive the corresponding kinetic equation for the dis- 
tribution function of this stochastic variable, i.e. the Fractional Fokker-Planck 
equation which has fractional space derivative instead of the usual Laplacian 
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(Sect.pT^). We show that the expression of the Fractional Fokker-Planck equa- 
tion is not unique (Sect.E^ and determine its scale invariance group, as well as 
its scaling solutions (ScctB). We also generalize the Einstein relation between 
the statistical exponents and the diffusion coefficient (Sect||). This helps us to 
clarify the pf^sical meaning of the exponents characterizing a stable distribu- 
tion. In Sect.g, we compare our approach to those followed by |3[ [l7| . 
Finally, we discuss (Sect.0) the possible applications of the theoretical results 
that we have obtained. 



2 Fractional Fokker— Planck equation 

2.1 Generalized Langevin equation 

We start with the Langevin-like equation for a stochastic quantity X{t): 

In the classical theory of a Brownian motion, X{t) is theJocation of Brow- 
nian particle under the influence of stochastic pulses Y{t) a. The statistical 
properties of this stochastic forcing will be specified below. We first need to 
derive an equation for the distribution function 

p{x,t)^{S[x-X{t)]) (3) 

where the brackets (...) denote statistical averaging over stochastic force reali- 
sations. Due to the fact that the Dirac function is the Fourier transform of the 
unity, we have: 

/°° rik 
—exp{-tk[x-Xm (4) 

When averaged, Eq.^ yields merely that the probability is the inverse Fourier 
transform of the characteristic function Zx{k,t) (see the Appendix A for an 
alternative derivation exploiting more directly this property): 

Zxik,t) ^ (expiikXit)])) (5) 
p{x,t) = F-^[Zxik,t)] (6) 



^ Following the approach of Einstein and Schmoluwski we neglect the inertial term for large 
time lags, and therefore consider the balance between viscous friction and forcing. 

■^These pulses correspond to a force divided by the friction coefficient, i.e. the inertial mass 
divided by the viscous relaxation time. 
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where F and F ^ denote respectively the Fourier-transform and its inverse: 



F[f] - f{k) = I dx exp{tkx)f{x) F-^[f] = f{x) = / ^ exp{^tkx)f{k) 



2tt 



On the other hand, Eq.0 can be integrated into: 



(7) 



dTY{T) 



(8) 



Since we can assume i without loss of generality that X{0) ~ 0, we obtain the 
following equation: 



dp 
dt 



F 



exp 



ik / dry(r) 



(9) 



Now, to make a further step, it is necessary to specify the statistical prop- 
erties of the stochastic source. We consider the particular example when 
the source is represented. as a sum of independent stochastic "pulses" acting at 
equally spaced times tj □: 



Yit)=J2Y,,AASit-t,) . (10) 

where to = 0, i^+i — tj — A (j = 0, 1, 2, ....) and the pulses Yj are independent 
stochastic variables having stable Levy distribution -P{Y,.a} for all j and which 
has the following characteristic function p| 



Zy^^Ak) = {exp{ikYj^A)) = expA ^^i-fk - D\k\'' 



(11) 



where a, /?, 7, D are real constants (0 < a < 2, — 1 < /3 < l,D > 0) and uj{k, a) 
is defined as: 



a 7^ 1 : uj[k, a) — tan—^; 



1: cjik,a)^-log\k\ (12) 



a and /3 classify the type of the stable distributions up to translations and 
dilatations: with given a and /?, 7 and D can vary without changing the type of 
a stable distribution. The parameter a characterizes the asymptotic behaviour 
of the stable distribution: 



p{x) 



00 



(13) 



■^Indeed, we are considering only the 'forward' Fokkcr-Planck equation. 
*See Sect.H for an alternative corresponding to a power-law distribution of the waiting 



times 
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hence, corresponds to the critical order of moments for their divergence: 

H>a: (x^) = oo, (14) 

For (additive) walks a is also related to the fractal dimension of the trail 
[ p!o| , whereas for the generator of the (multiplicative) universal multifractals it 
measures their multifractality The parameter f3 characterizes the degree 

of asymmetry of distribution function. Indeed, if /3 = 0, then negative and 
positive values of Yj^a occur with equal probabihties, while if /? = 1 or /3 = — 1 
(maximally asymmetric distributions) then, for < a < 1 and 7 = P{Yj,a} 
vanishes outside from [0, +00] or respectively from [—00, 0] u. We already men- 
tioned that maximal asymmetry is required for generators of universal multi- 
fractals; let us add that in this case the Laplace transform is more convenient 
than the Fourier transform. The nonzero value of P implies the existence of 
a primary direction of the stochastic pulses (that is, the direction to plus or 
minus infinity), and thus the existence of a drift for particles in this direction. 
For more details concerning the properties of stable laws see, e.g. The 
meaning of 7 and D will be discussed and clarified below. 

Now, using Eq.^ and the independence condition of the stochastic pulses 
Yj,A we get: 



= ]^(ea;p(ifcYj-,A)) = {exp{ikYj^A)) 
3=0 

(15) 

where n is a number of pulses corresponding to the present time t — nA. 
Therefore , with the help of the equation of the characteristic function of the 
pulses (Eq.|l^), we obtain the characteristic function Zxik,t) ( Eq.^ of the 
stable process: 



exp 



ik I dry^r) 




exp 



3=0 



'Y, 



3A 



Zx{k,t) 



exp 



ik / dTY{T) 





exp |i 



i-ik-D\k\°' ( 1 - i(3—uj(k,a) 



(16) 

The fact that this process has stationary independent increments (i.e. 
pulses 5^-, a) gives the possibility to get directly Eq.|l^ without using any dis- 
cretisation of Y{t) as previously done (Eq.p^). Such a derivation is presented 
in Appendix A. 

Now inserting this expression of Zx(k,t) into Eq.||, one obtains: 



dp 
di 



dk 
2tt 



i-fk - D\k\°' + iPDuj{k, a)k\k\°'-^]Zx{k, t)exp{~ikx) (17) 



'For a > 1, P{Yj ,\} decays faster than an exponential on the corresponding half axis. 
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For the sake of the simplicity of notations, we will consider in the following 
only the case a 7^ 1, or (3 ^ 0. Therefore, Eq.^ reduces to: 

TTQ; 

Lo{k,a) = Lo{a) = tan-^ (18) 

2.2 An expression of the Fractional Fokker— Planck Equa- 
tion 

One can see that in Eqjl^ the following type of integrals appears F~^{\k\°'Zx], 
which in fact correspond to fractional differentiations. Indeed, one may use 
Laplacian power for the Riesz's definition of a fractional differentiation since for 
any function f{x): 

- Afix) = F-\\k\'f{k)) (19) 
yields a rather straightforward extension: 

(-A)"/2/(x) = F-\\krf{k)) (20) 

Then, Eqjl^ yields: 



dp_^ dp ^ _^ 
dt dx 



{-Ar/'p + /3c.(a)^(-A)("-i)/2p 



(21) 



which for symmetric laws P — is a straightforward generalization of the clas- 
sical Fokker-Planck equation, by: 

A^-(-A)"/2 (22) 

This also points out that the scale parameter D of the Levy distribution 
corresponds to the diffusion coefficient of the Fractional Fokker Planck equation. 
On the other hand, the second term in the left hand side of Eq.|2l| has an obvious 
physical meaning. Independently on the value of a, it describes the convection 
of particles by the (constant) velocity 7. For a > 1, 7 corresponds furthermore 
to the mean value of the source {Y{t)), whereas it is no more the case for a < 1 
since the latter is no longer finite. In the latter case, the diffusion term has a a 
derivation order smaller or equal to the convection term. This confirms that the 
case a = 1 is indeed critical between two rather distinct regimes and it is more 
involved than other cases. Besides, it is worthwhile to note the role of the term 
(on the r.h.s.) related to asymmetry (f3 ^ 0). On the one hand, this term can 
be interpreted as an additional contribution to the convection due to existence 
of the preferred direction of the pulses related to {(3^0). On the other hand, 
such a flow is not proportional to p (as the convective flow does) but rather 
to (— A)("~^^/^p, which is rather typical for the diffusion flow. In some sense, 
due to this term the division of flows into convective and diffusion ones (as 
done in the standard Fokker-Planck equation) becomes rather questionable and 
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presumably no longer relevant for the Fractional Fokker-Planck equation. One 
may note that a somewhat similar weakening of this distinction occurs also in 
the classical Fokker-Plank for nonlinear systems On the other hand, it is 
easy to check that the Fractional Fokker-Planck equation is Galilean invariant, 
as is should be: the velocity of the moving framework just add to 7. 

2.3 The non uniqueness of the expression of the Fractional 
Fokker-Planck Equation 

One cannot expect to obtain a unique expression for the Fractional the Fokker- 
Planck equation, since there is not a unique generalization of the differentiation 
to a fractional order. Indeed, there exist various definitions of the fractional 



differentiation (see, e.g. |29 and references therein) which are not equivalent. 
This will be illustrated by two examples in the next section. The first one is 
related to the fact that there are 'signed' (fractional) differentiation and respec- 
tively 'unsigned' (fractional) differentiations, i.e. differentiations which are not 
invariant and respectively invariant with the mirror symmetry x —x. In the 
case of standard differentiation, the question of signs is fixed: 'signed' and 'un- 
signed' differentiations correspond merely to odd and respectively even orders 
of differentiation (hence the unique expression of the classical Fokker-Planck 
equation, which is of second order). This is no longer the case for fractional 
differentiations. 

The second example corresponds to the fact that fractional differentiations 
are in fact defined by integration, and therefore can depend on the bounds of 
integration. 

Nevertheless, we are convinced that the expression corresponding to Eqj2^ is 
at the same time the simplest one to derive and the one whose physical signifi- 
cance is the most straightforward. On the other hand, let us emphasize that the 
existence of distinct expressions for the Fractional Fokker-Planck equation does 
not question the uniqueness of its solution. Indeed, these distinct expressions are 
equivalent because their solution should correspond to the unique probability 
density function corresponding to a given Langevin-like equation (Eq.^). 

The non uniqueness could be rather understood in the following way: cor- 
responding to the distinct fractional differentiations (and their corresponding 
fractional integrations), there should be distinct ways of solving the Fractional 
Fokker-Planck equation in order to obtain its unique solution. 

2.4 Two alternative expressions of the Fractional Fokker- 
Planck Equation 

Contrary to the unsigned fractional power of a Laplacian Eqj2^, let us consider 
for instance the following 'signed' fractional differentiation: 
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^f{x)^F-'[i-^krm]. (23) 
With the help of (i) the identity {9{k) being the Heaviside function): 

|fc|" = fc"[6i(fc) + (-l)"6'(-fc)] (24) 
and of (ii) the inverse Fourier transform of the Heaviside function: 

^-M^W]4^(.) + ^ (25) 

as well as of (iii) the property that a Fourier transform of a product corresponds 
to the convolution of the Fourier transforms, one derives from Eq.|l^ an another 
form of the Fractional Fokker-Planck equation. 

9p dp „/ Tra . na TTa\ d°'p . , . ira 9" /" dx' p{x' ,t) 



dp dp f 7ra ira ■Ka\ d°'p , , -na 9" f 

m+^d^ = r^^- + Ps^n-tan-) —^D{l-(})s^n-— J 

n _ (26) 

Indeed, with the help of the following determinationaa (— i)" = e * 2 , (—1)" = 
g-^oTT^ Eq.|| yields: 

{~ik)°'[0{k)e''^ + e{-k)e-''^] (27) 



TT X — X' 



and with the help of Eqs.^,^ 27, it is rather straightforward to derive Eq.p6|. 

However, Eq.|2^ is already rather involved in the case /? = 0, whereas this 
case is obvious for the equivalent Eq.pl|: 

dp dp nad°'p na d"' ^ dx'p{x',t) 

the last term of the r.h.s. of Eq.^ is rather complex, whereas indispensable. 
Indeed, there is a need of signed second term to counterbalance the first signed 
term of Eq. in order that the r.h.s. of Eq|2| will correspond to an unsigned 
differentiation (the fractional power of the Laplacian in Eq.^l]). Both terms 
correspond to the signed fractional differentiation of order a but whereas it is 
applied to p in the former term, it is applied to an integration of a zero order 
of p in the latter term. This zero order integration corresponds to the effective 
interaction of particles having a scaling law inversely proportional to the distance 
between them. An analogy with the interaction between dislocation lines 
can be mentioned. It is plausible that the collective effect corresponding to 
this the effective interaction of particles could be responsible of the large jumps 
which are so important in Levy motions. 



'One may note that the existence of 



her determinations confirms the non uniqueness 
of the fractional derivative defined in eq.|24. Furthermore, taking another determination will 
merely modify some prefactors in r.h.s. of Sq.p^ 
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An other expression of the Fractional Fokker-Planck equation can be also 
obtained with the help of the Riemann-Liouville derivatives. The ^-th order 
Riemann-Liouville derivatives on the real axis are defined as 



^ ^ ' T{l-fi)dxJ_^ {x-x'Y ^ ' T{l-fi)dxJ, {t-x'Y 

(29) 

where , D^^ are respectively the left-side and the right-side derivatives of 
fractional order /i (0 < /i < 1) and T is the Euler's gamma-function. Appendix 
B gives a derivation of the corresponding expression of the Fractional Fokker- 
Planck equation which is: 

9i 9a; ^ ^ + p - f/Jwiaj— Dl p (30) 

3 Scaling properties of the Fractional Fokker- 
Plank equation 

Now let us return to one of the equivalent expressions (Eqs.^ |2^, 30) of 
the Fractional Fokkcr-Plank equation and consider its scaling group properties 
which defines fundamental properties of its solutions. Without loss of generality 
we can assume 7 = 0. The scale transformation group can be written in the 
following manner: 



t^Xt', x = X^x', D^X'^D' (31) 
p{x, t; a, 13, D) = X^p{{x', t'; a, /3, D') (32) 

Here Xi '^i S are yet unknown exponents of the scale transformations which 



wn e xp i 

should leave invariant the Eqs.|l], |30[ Onc_may note that the subgroup not 



dealing with D transformations was used by [g^ in order to obtain the fractional 
order of differentiation of the Fractional Fokker-Planck equation corresponding 
to a symmetric stable processes (see Sect.^. Due to the further normalization 
condition for the distribution function, one obtains the following 2-parameters 
scale transformation group: 

t = Xt', x = X^x', D = X'^^-^D' . (33) 

p{X''x,Xt;a,f3,X°"'-^D) ^ X~^p{x,t;a,(3,D) . (34) 

X, X being the arbitrary group parameters. 

If furthermore the initial condition p(x, 0) is invariant under the scaling group 
(Eq,33), then the corresponding solution of Eqs.|^, ^ remains invariant 
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under the action of this group for any other time. The simplest example of an 
invariant initial condition is p{x, 0) = S{x), which is of fundamental importance 
since it corresponds to the Green functions of Eqs.^, ^ 

Let us analyze the general properties of the invariant solutions in a similar 
way to renormalization group approach |3|] (analogous consideration was 
used in a more complex variant when calculating the spectrum of a compressible 
fluid iQ). Due to the fact that the scaling invariant solutions should satisfy 
the identity (Eq.^) for any value of the arbitrary parameters A, x, they could 
depend only on products of variables which are independent of them. Therefore, 
due to the relationships ax — k — 1 = and S + k — 0, the scaling solutions are: 

where $(.) — p(l, .), 5* = p{., 1) are arbitrary functions which are determined 
by the initial conditions. Therefore, Eq.^ represents the general form of the 
invariant solutions of the kinetic equation Eqs.|l], p6| , ^ 

Eq. ^ can be obtained by first differentiating Eqj34| with respect to A, and 
then, to X ^nd setting A = 1,X = — This yields the following system of 
equations: 

dp xdp 1 9P , r, cqR^ 

at a ax a at oD 

which are linear and therefore can be solved by the method of characteristics 
and their solution indeed correspond to Eq.Ka. 



4 Some particular solutions 

4.1 Explicit solutions 

With the exception of the three following cases, there is no way to obtain an 
explicit expression of the solutions, with the initial condition f(x,0) — 6{x), of 
the Fractional Fokker-Plank equation (Eqs.|l], ^ in a closed form with the 
help of elementary functions: 
1) a = 2 (/3 = 0): 

It corresponds to Gaussian distribution of the stochastic forcing and to the 
classical Fokker-Planck equation: 



dp ^ 9p ^ ^ 
dt dx 

which solution is the normal distribution: 



d'^p 
dx"^ 



(37) 



p{x,t) 



1 



exp 



- It? 
ADt 



(38) 
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2) a = l,/3 = 0: 

It corresponds to a forcing having a Cauchy distribution and to the following 
Fractional Fokker-Planck equation: 



dp_^ dp ^ ^d_ 
dt dx dx 



dx' p(x' , t) 



J TT X — X 

The solution of Eq.|39| with the initial condition f{x, 0) = 5{x) is: 

Dt 1 



p{x,t) 



TT {X-Jt)^+DH^' 



(39) 



(40) 



3) a = 1/2, /5 = 1. 

Then the expression displayed in Eqj2^ of the Fractional Fokker-Planck equation 
takes then the form: 



dp^ dp ^ d^/^p 
dt dx dx^l"^ 
The solution of Eq.^ with the initial condition /(x, 0) = &(x) is: 



p(x, t) = 9{x — 7<) 



Dt 



y/27r{x - -itf/- 



--exp 



DH 



2j-2 



2{x--tt) 



It is easy to check that all the explicit stable distributions (Eqs. 
belong to the class of scale invariant solutions. 



(41) 

(42) 
§,|2D 



5 Generalization of Einstein relation and anoma- 
lous diffusion coefficient 

The scaling analysis of the moments of the distribution function will lead to the 
generalization of the Einstein relation for the anomalous diffusion. However, 
there is an important difference, since moments of order larger than a < 2 will 
diverge and in particular: (x^) = oo. On the other hand, the motion remains 
mono-fractal, since all the moments (a;^) (0 < /i < a) will have the same scaling 
law. 

The statistical moments of the distribution function of particles initially 
concentrated at the origin (p(x, 0) = 5{x)) correspond to : 



dxx^ 



dk 
2^ 



exp 



-ikx — D\k\"t [ a — iP-rritO'n--^ 
\k\ 2 



(43) 



In agreement with the scaling properties obtained in Sect.^), we have: 

r^(t) = (a;^)!/^ - {aDty/"C{a, (3, p) (44) 
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which is obtained by renormahzing x and k by (_Di)^/", i.e. by considering the 
following variables: 

x,=k{aDtY/-^X2 = j-^ (45) 

which yield from Eq.^ the following prefactor (depending neither on time nor 
on LI) : 



C(a, Am) = I r ^...^ r f-i-^ - (l - ^Ps,nix,)tanC-^) 

U-oo J-ao 27r [ a \ 2 

(46) 

It follows from Eqjij that the scaling of rp(t) in respect to D and t is 
universal and does not depend on the order /i of the moment considered (/i < a). 
Therefore the Einstein relation can be formulated in terms of any of the finite 
moments. Indeed, only the the numerical prefactor C{a,l3,^) (Eq.^) depends 
on fi, but neither on time nor on D. 

The Gaussian case yields the classical Einstein formula : 

r2(t) - i2Dty/^ (47) 

Not surprisingly, Eq.^ confirms the fact, already pointed out on Eq.|l], that 
D does corresponds to a (generalized) diffusion coefficient. 

On the other hand, let us confirm that the scaling behaviour (Eq.^) is 
independent of the initial distribution of the particles. This is due to the fact 
that the distribution with initial condition f(x, 0) = S{x) plays the role of the 
Green function for the distributions with other conditions, and therefore imposes 
its scaling on time and on D. This is easily confirmed with the help of Eq.^ 
which gives the general expression of the scaling probability densities: 

with: 

f roc ^ 



6 Comparison with other approaches 

As mentioned in Secul, particular cases of the Fractional Fokker- Planck equa- 
tion were obtained ]12| , 13| , Q , and on the other hand a Langevin-type equation 
has been obtained |17|] for the nonlinear Fokker-Planck equation |l5l uM whose 



solutions maximize the generalized q-entropy introduced by |18 and exhibit a 
Levy-like anomalous scaling. 
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A Fractional Fokker-Planck equation was obtained by ||l2|, |l^ in the frame- 
work of the continuous time random walks (CTRW's) model of anomalous dif- 
fusion |35j. However, this method does not involve directly a Levy process, but 
a walk sharing some common behaviour of the latter, without being equivalent 
to it. Indeed, the distribution of steps, which corresponds to the probability 
distribution of the pulses in the Langevin equation, is considered as a pure 
power-law. This corresponds only to the asymptotic behaviour of the Levy dis- 
tribution, i.e. its tails, and therefore takes into account only one of the Levy law 
parameters. This nevertheless allows to establish scaling relations, therefore to 
determine the fractional order of differentiation (see however a remark below), 
but not to determine a precise expression of this fractional differentiation. This 
is already the case for its coefficient, i.e. the (fractional) diffusion coefficient 
of the Fokker-Planck equation, since scaling reasoning does not yield a relation 
with the scale parameter. Second, there is no simple way to deal with the skew- 
ness parameter (3 when considering the probability distribution. Therefore, the 
corresponding non trivial term in the Fokker-Planck (Eqs.pl], p6| , 30) was not 
obtained by | [T^ , [T^ . On the contrary, in our approach the four parameters, 
which determine all the statistics of the pulses, are all taken into account in 
an exact and rather straightforward manner with the help of the characteristic 
function. On the other hand, jl2, 13 show an easy generalization to both tem- 
poral and spatial memories is obtained by introduced a a second (generalized) 
Langevin equation for the waiting times, which introduces a Fractional Fokker- 
Planck equation for the latter. However, let us point out that the usual scaling 
reasoning does not apply in a straightforward manner due to the divergence 
of moments associated to the power-law of the Levy probability distribution. 
Indeed, one cannot consider the scaling of the variance of the distance (r^(s)) 
traveled by a particle after s steps, because it is a divergent statistical moment, 
i.e. equals to infinity, as soon as a 7^ 2. One has to consider moments of order 
fj, < a, which are not only finite but are furthermore monoscaling (Sect. 
This last property means that the scaling of the fi — th root of moment of order 
fi is independent of this order, and therefore explains why the (mono-) scaling 
reasoning works. 

In a recent article |lj] a different form of fractional Fokker-Planck was in- 
troduced with the help of a phenomenological and interesting modification of 
the classical Fick law into a fractional Fick law, which is discussed in details. 
The question of the existence of a corresponding Langevin-like equation remains 
open. Some integral convergence problem imposes that the power-law exponent 
of the probability distribution tails belongs to ]1, 2]. This exponent would corre- 
spond to the Levy stability index a the solution of this fractional Fokker-Planck 
was a Levy motion. However, this is not exactly the case, although it seems at 
first glance rather similar to it. Indeed the corresponding characteristic function 
(see Eqs. 14-15 of Ref. involves k" instead of (ifc)" for the characteristic func- 
tion of a Levy motion (Eq. |lT| ). The difference is more obvious when considering 
the asymmetric extension, either on the proposed Fokker-Planck equation (Eq. 
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17 of Rcf. iQ or on the characteristic function (Eq. 19 of Ref . ) , since it does 
not include the non trivial asymmetric term that we put in evidence (Eqs.^, 
p6| , |30| ). Nevertheless, it corresponds to an interesting variant of asymmetric 
diffusion, which solutions could be rather close to stable Levy distributions. 

An essentially different generalization of the Fokker-Planck equation was ob- 
tained by |l^, |l^ for anomalous diffusion. Indeed, |l^ showed that the solutions 
of the nonlinear Fokker-Planck equation introduced by |l^, |l^, maximize the 
generalized q-entropy ||l^ , and correspond to a well defined Levy-like anomalous 
diffusion, but with a finite variance and non zero correlation. Both properties, 
which seem relevant and desirable for many applications, in particular for the 
transport in porous media, are not satisfied by a Levy motion. 

Furthermore, demonstrated that the corresponding Langevin-like equa- 
tion has the particularity that the random forces are modulated by a given power 
of the probability distribution. This corresponds to a macroscopic feedback to 
the microscopic kinetics, which is absent in our Langevin equation. 

In comparison with these different works, we followed a rather distinct ap- 
proach since we started with a Langevin-like equation with random forces which 
are exact stable Levy processes, which can be symmetric as well as asymmetric, 
and with no limitation on the possible values of the Levy index a. The partic- 
ular case corresponding the symmetric stable processes was previously inferred 
by [|l2[|l|,[l4|. However, we showed that in the more general case of asymmetric 
stable processes, a new non-trivial term appears which has a rather intermedi- 
ate role between diffusion and convection (see Sect. 2.2). Firrthermore, the use 
of the characteristic function allowed us to obtain a generalization of Einstein's 
formula. We also clarify the fact that different expressions of the same Frac- 
tional Fokker-Planck equation are obtained, depending on the type of fractional 
differentiation which is used. 

On the other hand, the conjecture issued by ||l6| that 'further unification 
can be possibly achieved by considering the generic case of a nonlinear Fokker- 
Planck -like equation with fractional derivatives' should be closely examined, as 
well as the fractional time evolution suggested by Fogedby, Compte. 



7 Examples of applications 

7.1 Generalisation of the Holtzmark distribution 

In the introduction of this paper, we recalled that the gravitational force result- 
ing from randomly and homogeneously distributed point masses which acts on a 
given test point mass has 0] has a stable symmetrical law with a — 3/2. One 
may note that a similar problem arises with charged particles and electrostatic 
forces. However, the distribution of masses in the Universe is rather inhomoge- 
neously distributed, e.g. on a fractal set of fractal dimension Q Dp sa 1.2. (see 
for discussion and a multifractal analysis). Let us extend the original result 
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of Holtzmark to the case of particles distributed on a (possibly fractal) space 
of dimension d and which mutually interact according to a scaling law L 
being the distance between two particles, e.g. the collective force f acting on 
the randomly chosen particle of unit mass located at x, which results from the 
distribution of masses rrifc at points x^*^^ , has the following type □ : 

fw-E-^- i/_''I(;)i"+i (49) 

Due to the linearity of Eq.^, the superposition ([m+m']) of two independent 
mass distributions {[m] and [m']) yields a force having the same probability 
distribution as the one of the sum of the two forces resulting from each of mass 
distribution, i.e.: 

i{[m]) + {{[m']) = {{[m + m']) (50) 

on the other hand, the fact that the mass is concentrated on a fractal set of 
dimension Dp, it should scale in the following manner with the (space) scale 
resolution A, i.e. the ration A = j of the outer scale L over the inner scale I of 
the fractal set, in particular in the limit A — > cxd : 

mx = miX-^'^ (51) 
which implies with the help of Eq.^ the following scaling for the forces: 

h[mi] 4 A^ViK] ' (^^) ^ "/iK] (52) 

which together with Eq.|50| demonstrates |Q that f\[m\\ has a (symmetric) Levy 
stable distribution, with a Levy index a — Dp/r. 

With this Levy index value for the random forces, and neglecting their inter- 
relations (which will be studied elsewhere), we may define the random velocity 
of the test mass as defined by Eqj| and its probability distribution by Eq.^ 



7.2 The anomalous diffusion of a passive scalar by a two- 
dimensional turbulence 

Let us consider the velocity field Wi(x) resulting from point-like vortices: 

where k is the intensity of the vortices, eij\ij — 1,2 is the fundamental anti- 
symmetric tensor, x*^ is the location of the fc*/i vortex. Following |pa we assume 



^by 'type', we mean that most of the algebraic details of the following equation are irrele- 
vant, only its scaling properties are relevant. 
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that the vortices are distributed on a fractal set of fractal dimension Dp. We 
are therefore in the situation of the generalization of Holtzmark distribution and 
indeed Eq. |5^ is of the type of Eq. with r — 1. Therefore, v[(k)] has a Levy 
probability distribution with a = Dp. one may note that this result can be ob- 
tained by using straightforward calculations of the characteristic function in the 
manner analogous to The main distinction is that the fractal distribution 
of the vortices had be taken into account. 

Therefore, the diffusion of a passive scalar in 2D turbulence created by 
a fractal set of point-like vortices is defined by Eq.|l| with ^ — [3 — Q and 
a = Dp. It points out the interest of using Fractional Fokker-Planck equation 
for the analysis of the diffusion processes of particles in turbulent media. 

7.3 Multifractal modeling 

However, as noted in the introduction, the most appealing area of application 
of the Fractional Fokker-Planck equation could be for 3D turbulence. Indeed it 
should play a key role for the definition of the generators of dynamic universal 
multifractals js^. Let us first recall some basic features of static universal 
multifractals, i.e. defined only on space. The corresponding field, e.g. the 
flux of energy Fx at higher and higher resolution A = j, should respect the 
multiplicative property of the scale ratio, i.e.: 

FA^x-y = Fx-Tx{Fx') (54) 

where Tx is a scale contraction operator of ratio A, which in the simplest case 
is the isotropic self-similar contraction (Tx(x) = A~^x). Therefore Fa might be 
defined with the help of the generator F of this group , more precisely speak- 
ing by the exponentiation of the latter which satisfies the following additive 
property: 

La^a.a' =rA + rA(rv) (55) 

In order to satisfy the multiscaling power law: 

VAe (1,A) : {F^)^X^'^'i^ (56) 
the generator should have a logarithmic scale divergence: 

La ~ log A (57) 

this latter condition is obtained by a convolution of a given the Green's function 
g over a white-noise 7a (called the 'sub-generator'): 

rA=.9*7A (58) 

In the case of universal multifractals jioj the sub-generator is an extremely 
asymmetric and centered Levy stable with a Levy index a and the condition of 
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logarithmic of divergence (for a Z?— dimensional isotropic process) corresponds 
to: 

5"fe)«Ur"-^«;i?H = - (59) 

a 

However, in order to take into account the causality for time-space processes 
[ p4| , it is rather more interesting to consider a differentiation operator, i.e. to 
consider g^^ rather than g, i.e.: 

g-\x,t)*rAix,t) ^-fAix,t) (60) 

Furthermore, in order to take into account the difference of scaling in space 
and time, the time and respectively space orders of differentiation should be 
different. Therefore, the following type of differential equation were considered: 

ff-^=| + (-A)i-^' (61) 

where the 'elliptical dimension of the space' D^i = D + 1 — Ht is the effective 
space-time dimension, D is the dimension of the space cut and Ht corresponds 
to the deviation of the time scaling in comparison to the time scaling. The 
Fractional Fokker-Plank equation (Eq.^lj) suggests that the following fractional 
operators could be as well considered: 

= I + (-A)i-^' +/3c.(a)A(_A)(i-«.-i)/2 (52) 
ot ax 

and for any value of P the evolution of the generator has a microscopic inter- 
pretation with the help of the corresponding Langevin equation, i.e. it points 
out that Eq.^ corresponds to a (generalized) path integral. 



8 Conclusions 

The original results obtained in this paper are the following: 

1. the Fractional Fokker-Planck equation, i.e. the kinetic equation describing 
anomalous diffusion in response to a stochastic forcing having a Levy 
stable distribution, which can be symmetric, as well as asymmetric, 

2. the a physical interpretation of all the parameters of the Levy stable distri- 
butions, due to a precise determination of the coefficients of the Fractional 
Fokker-Planck equation, 

3. the scale transformation group of the Fractional Fokker-Planck equation, 
as well as corresponding scaling solutions, 

4. the universal dependence of the distribution function moments on the 
diffusion coefficient and time. 
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5. some preliminary examples of applications, including a generalisation of 
the Holtzmark distribution, two-dimensional diffusion of a passive scalar 
and multifractal modeling of intermittent fields. 

We compare these results with particular cases obtained by |l|, |3[ |l|, 
as well as their relations to a nonlinear Fokker-Planck equation introduced by 
^ whose solutions exhibit an anomalous Levy- like diffusion ||l^, [l7|] . 

In summary, we believe that the kinetic equation obtained will be useful for 
studying various physical systems with non-Gaussian statistics. 
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A Another Way to Derive the Fractional Fokker- 
Plank Equation. 

Here we present another approach to the derivation of Eq.n6l Let be the prob- 
ability distribution function p{x, t) of x{t) which obeys Eq!g and having first 
characteristic bmction Zx{k,t): 

Pix,t) = F-'[Zx{k,t)] (A.l) 

Due to the linearity of the Fourier transform and the fact that the one 
considered applies only in space, not in time: 

dp^^_,(dZxik,tn ^^^^^ 



dt V dt 

If the pulses Ydt{t)' s for infinitesimal time lag dt are independent and iden- 
tically distributed variables for any arbitrary time t, and have a second charac- 
teristic function dtKyik): 

ZYdt{k,t) = exp[idtKY{k)] (A.3) 
then X{t) has independent increments [2^ and: 
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Zx{k,t) = exp[itKY{k)] (A.4) 

The demonstration of Eq. A.4 is rather straightforward, but can be also ob- 
tained with the help of the time discretisation which we used in order to obtain 
Eq.^ 

Inserting Eq.(A.4) into Eq.(A.2) we get 

= F-'[^Kyik)Zxik,t)] (A.5) 
The particular case of Levy stable pulses Y{t) corresponds to 



Kvik) = -fk + iD\k\°' 



k 

1 - i0jj^^^^{k,a) 



(A.6) 



therefore, we immediately get Eq.|T^, and we only need to interpret |A;|" and 
in the physical space, as done in Sect. 2. 



B Fractional Fokker— Planck Equation in Terms 
of Riemann— Liouville Derivatives. 

The /i-th order Riemann-Liouvillc derivatives on the infinite axis are defined as 

where D'^,D^ are respectively the left-side and the right-side derivatives of 
fractional order /i (0 < /i < l),r) is the Euler's gamma-function. 

The Fourier-transforms of the fractional derivatives (B.l), (B.2) are the 
following: 

i^(D^/) = (T^krfik) (B.3) 
where /(fc) is a Fourier transform of /(x), and we have: 

(Tik)" = \k\''exp (^T^signk^ (B.4) 
where jj, is real and which yields: 
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therefore 



{^.^kY{^+ikYfk ^F-^ fc2^/fe IfcpA'/fe (B.5) 



(B.6) 



Eq. B.7 establishes the equivalence between Eq.pT| and Eq.|30|. 



(B.7) 
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